A thermodynamic characterization of future singularities? 
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Abstract 

In this Letter we consider three future singularities in different Friedmann-Lemaitre- Robertson- 
Walker scenarios and show that the universe departs more and more from thermodynamic 
equilibrium as the corresponding singularity is approached. Though not proven in general, this 

q h feature may characterize future singularities of homogeneous and isotropic cosmologies. 
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I. INTRODUCTION 



Cosmological singularities arise as displeasing features in mathematical models of the uni- 
verse [l[ S|; world lines terminate and/or physical quantities diverge. It is usually argued 
that this results because of the exceedingly simplicity of the models in question, implying 
that when we reach a much better understanding of the physical processes that take place 
under the most extreme conditions we will be able to design realistic models free of singu- 
larities. 



Perhaps the most worrying of all is the Big Bang singularity which persists even if the stan- 
dard cosmological model is corrected with the addition of an era of inflationary expansion 
just before the radiation dominated epoch. This explains the interest raised by proposals of 
universes, such as bouncing {3] and emergent with no beginning at all. 



Singularities are heralded, among other things, by the growing without bound of key physical 
quantities, such as energy densities and pressures. Here we focus on future cosmic singular- 
ities in Friedmann-Lemaitre-Robertson- Walker (FLRW) universes and ask ourselves if they 
may be characterized by some thermodynamic distinctive feature. Our provisory answer is 
in the affirmative. We reach this tentative conclusion after considering three singularities, 
linked to different cosmological models. As it turns out in the three cases the total entropy of 
the systems does not tend to a maximum as the singularity is approached. That is to say, the 
systems do not move towards thermodynamical equilibrium but on the contrary: the closer 
they get to the singularity, the further away from the said equilibrium the systems depart. 
By "system" we mean the apparent cosmic horizon plus the matter and fields enclosed by it. 



At this point it seems suitable to recall that physical systems tend spontaneously to some 
equilibrium state compatible with the constraints imposed on them. This summarizes the 
empirical basis of the second law of thermodynamics. Put briefly, the latter establishes that 
isolated, macroscopic systems, evolve to the maximum entropy state consistent with their 
constraints [5]. As a consequence their entropy, S, cannot decrease, i.e., S' > 0, where the 
prime means derivative with respect to the relevant, appropriate variable. Further, S has to 
be a convex function of the said variable, S" < 0, at least at the last phase of the evolution. 
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The apparent horizon in FLRW spacetimes is defined as the marginally trapped surface with 
vanishing expansion of radius js] 



1/VH 2 + ka~ 2 , (1) 



where a and k are the scale factor of the metric and the spatial curvature index, respectively, 
and it is widely known to have an entropy which, leaving aside possible quantum corrections, 
is proportional to its area 

S h cc A = 4vrf^ , (2) 

which agrees nicely well with the holographic entropy derived from considerations of the 
foamy structure of spacetime [7]. Besides, this horizon appears to be the appropriate ther- 
modynamic boundary jg]. 

In its turn, the entropy of the fluid enclosed by the horizon is related to its energy density 



and pressure by Gibbs' equation [5], namely, 

T dSf = d (p ^f^j +pd f^-f\\ , (3) 
where T stands for the fluid's temperature. 

As said above, for the second law to be satisfied the inequality dSh + dSf > must hold at 
all times, and d 2 (Sh + Sf) < at least at the last stage of evolution. 

Application of this broad idea to cosmic scenarios leads to a variety of interesting results. 
Among others, dark energy (or some other agent of late acceleration) appears thermody- 
namically motivated: both in the case of Einstein gravity [9|] and in modified gravity [10 ]. 
In particular, ever expanding universes dominated either by radiation or pressureless mat- 
ter cannot approach thermodynamic equilibrium at late times. This is also true for those 



phantom dominated universes whose equation of state parameter, w = p/p, is a constant js]. 



The target of this Letter is to study whether the universe gets closer and closer to thermo- 
dynamic equilibrium as it approaches a future singularity. We assume Einstein gravity and 
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provide some general relations in section [TTJ Then we consider three specific FLRW scenarios: 
the big crunch singularity (section UTB . a sudden singularity (section ITVl) . and a "little rip" 
singularity (section [V}. Discussion and final comments are presented in section IVT1 



II. GENERAL RELATIONS 



The field equations for a spatially homogeneous and isotropic universe are the Friedmann 
equation 



3 ( if 2 + 4 I = 4 = 87rGp, (fc = 0,±l), 



and 



H =-4ttG (p + p) + ^ = -H 2 --^ [1 + 3- 

cr 2r\\ p 



P 



Quite generally we find 



a V P 



where the prime denotes the derivative with respect to the scale factor a, and 



A" = 36tt4 + - 
a 2 \ p 



(4) 
(5) 
(6) 

(7) 



Via the proportionality Sh oc A (cf. <^), the relations © and (j7]) will allow us to obtain 
information about S' h and S'^, respectively, for various choices oip/p and p'/p'. 

To determine the derivatives of Sf we must first discern the temperature evolution of the 
fluid. From Euler's relation nT s = p + p, where n and s are the number density of parti- 
cles in a comoving volume and the entropy per particle, respectively, and the conservation 
equations p' = —3(p + p)/a and n' = —3n/a, we find 

s'=p'-(p + p)y ■ ( 8 ) 
Taking into account the perfect fluid condition s' — 0, the temperature behavior is governed 
by 

V v' 

= ■ (9) 

T p + p K ' 

Straightforwardly one obtains the desired expressions, 



S'=2n 



?aP + P 



and 



r 3 A P + P 
T 



Q/l 



2tt 



12 



P 



1 + 



3 p 
2p 



1 + 3 



P 



p' 



i + 



p 



+ 



1 + 3 



P 



(10) 



(11) 



Sections III to V apply the above set of formulas to three different cosmological scenarios 
that harbor a future singularity. 

III. BIG CRUNCH SINGULARITY 

Here we focus on the radiation-dominated, spatially closed (k = +1) FLRW universe and 
explore the entropy behavior of the apparent horizon and radiation enclosed by it as the 
big crunch draws close. Under such circumstances the temperature gets so high that mat- 
ter becomes extremely relativistic and behaves as radiation, hence p = p/3 and pa 4 = const. 

The scale factor and Hubble function are given by 

*)=^l-(l-^)" * S f = -^, d2) 

where t s is the time at which a = and C 2 = ^p- pa 4 . The horizon area is A = ^ a 4 . 
Either by direct calculation or as special cases from <^ and ((7j) with p = p/3 it follows that 

^' = ^a 3 >0, and ^4" = ff a 2 > , (13) 

respectively. Thus, the graph of A increases and is concave for all values of the scale factor. 



Equation (jTUJ) leads to S'j oc a 2 > and equation ffTTj) implies oc a > . Thus, the 
total entropy (that of the apparent horizon, Sh, plus that of the radiation fluid), fulfills the 
generalized second law, S' h + S'j > 0. However, because S'^ + S'j > the system does not 
tend to thermodynamic equilibrium as the big crunch is approached. 
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IV. SUDDEN SINGULARITY SCENARIO 



In this scenario the singularity occurs at finite time (say t s ) and is characterized by the 
divergence of the acceleration and pressure while the energy density, scale factor and Hub- 
ble expansion rate remain finite. As the singularity is approached, the dominant energy 
condition (| p \< p) is violated but all other energy conditions are respected. 



et us consider the scale factor of the spatially flat (k = 0) FLRW metric as introduced in 
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121 ]. namely: 



tg / \ te. 



<*) = r - 1) + 1 - 1 - - , (14) 



where a s = a(t = t s ) > 1, and a and /3 are constant parameters lying in the range (0, 1] and 
(1, 2), respectively. Obviously, this expression holds for < t < t s . 



Since 

^ = t.^{tT + t. Hf'* - (15) 

the Hubble function H = a/a never becomes negative. 



A subsequent derivation yields 



m = j 2 (a - 1) (a s - 1) fpj a 2 - (l-f) 2 ■ (16) 



s 



In the limit t — >■ t s the first term, both in a(t) and a(t), dominates over the second one 
while the latter dominates in the expression for the acceleration which becomes negative 
and diverges. As a — > a s , H — > H s and p — > p s > where a s , H s and p s are all finite but 
p s — > oo via the field equation 3'd/a = —4ttG(p + 3p). 

In view of the above, in the said limit (t — > t s ) we can write 

a(a) = £ (a s - 1)1 , H(a) = f (a s - 1)^ , (17) 

t s \a s -lj t s a 

where we have eliminated the cosmic time in favor of the scale factor. 
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The area of the apparent horizon is given by 

A-tolT***, (|)V !)- %_£..„;,, ■ (18) 

where the second equality holds to leading order only. From ([6]) one has 

A! > (19) 

always. This corresponds to H'(a) < 0, which can also be checked explicitly from the full 
expression of H(t). Since p' > (the pressure diverges) and p' < 0, both terms in the square 
bracket on the right-hand side of equation (J7J) are positive, consequently, 

A" > (20) 

as well. Both A' and A" diverge upon approaching the singularity. 

Before concluding that in this sudden singularity scenario the universe departs from ther- 
modynamic equilibrium we must, as before, examine the thermodynamic behavior of the 
gravity source. 

From (11 01) and fill I) it follows immediately that with p > 0, p' > and p' < one has both 
S'j > and S'j > 0. The inequality S'j > together with A' > implies that the generalized 
second law, S' h + S'f > 0, is fulfilled in this scenario. On the other hand, the inequality 
S'h + S'j > means that the total entropy is a concave function of the scale factor and the 
universe does not tend to thermodynamic equilibrium. 



V. LITTLE RIP SCENARIO 

The expression "little rip" was coined as a contraposition to "big rip" [13]. In this scenario 
the ratio p/p < —1 but it increases and approaches —1 as time goes on. Although sooner 



or later all bound structures rip apart, at variance with the usual big rip scenario [131 ]. 
neither the energy density nor the scale factor diverge at a finite time. There is a future 
singularity but, because the expansion rate approaches de Sitter, it is pushed to t — > oo [l^ ]. 
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The equation of state of the gravity source is p = —p — Ap 1 / 2 , with A > 0. This alongside 
the conservation equation p' = —3(p + p)/a produces 

•L.4. / (/ \ | 1 

P = Po 



Inl- 



+ 1 



(21) 



. 2 \/Po \a>o, 

(cf. Eq. (12) in [14] ) . Notice that p augments with expansion, a typical feature of phantom 
dark energy, but logarithmically only and diverges just when a — » oo. The zero subscript 
denotes some reference time, we conveniently take it as the time at which the dark energy 
overwhelms all other components (matter, radiation, etc) to the point that their dynamical 
influence can be safely ignored. 



Because the FLRW metric is spatially flat, we have A = 4it/H 2 oc 1/p. From fl6]) one verifies 
A' < and ((7|) provides us with A" > 0. 



We next study the entropy evolution of the fluid. The equation of state p=-p- Ap 1 ' 2 in 
( ITUl) yields S 1 * > 0. Since A' and S'j have opposite signs, a closer look at the behavior of 
these quantities as they approach the singularity is necessary. Explicitly, we have 

.4 



A' 



-12tt 



and 



aH 3 T 



P 



1/2 



(22) 



(23) 



To proceed, information about the dependence of the temperature on a is required. The 
general law © yields 



T c ' 



.4 



2pV2 



(24) 



This integrates to 



T 
T 



3A f a . 
In — ) + 1 

2 \/po V a o 



1/2 



Consequently, 



31/2 2 + 4 
S' f = 2~4 ' n/ n 



34. 

/2 



(25) 



(26) 



T a^H 3 

In approaching the singularity, p diverges and the term oc p~ 1 / 2 in the second numerator of 
( 126]) can be neglected. Since H oc p 1 / 2 , one has S'j oc l/(a 4 p 3 / 2 ) in the limit of large values 
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of a. In the same limit A' in ( I22p behaves as A' oc — l/(ap 3 / 2 ). Therefore, as a — > oo the 
ratio | A'/S'f | diverges as a 3 whence S' h + S'j- < in the long run; i.e., the generalized 
second law is violated as the little rip singularity gets closer and closer. 



Despite we already know that in this scenario no thermodynamic equilibrium is approached 
when nearing the singularity, we shall determine whether the function Sh + Sf is convex or 
concave in that limit. From (J7J) we find that A" oc l/(p 3//2 a 2 ) for a — > oo while (iTTj) yields 
Sf oc — l/(p 3//2 a 5 ) in the far-future limit. The ratio | A" / S'j | diverges with the same power, 
a 3 , as the corresponding ratio of the first derivatives, i.e., S% + S'l > as a — > oo. 



VI. CONCLUDING REMARKS 



Macroscopic systems moving by themselves away from thermodynamic equilibrium is some- 
thing far removed from daily experience. This is enshrined in the second law of thermo- 
dynamics that introduces the entropy function and dictates its overall evolution. In this 
Letter we have studied the behavior of three FLRW universes as they draw close to a future 
singularity. As long as the entropy concept is related to the apparent cosmic horizon, none 
of the three approaches equilibrium. 



Against this some comments may be raised: (z) The proposed universes look too academic; 
we should not wonder that unrealistic systems do not comply with thermodynamics, (ii) 
The Universe is a very particular and unique system; why should it obey the thermody- 
namical laws? (iii) Everyone expects the breakdown of physical laws at singularities. Then, 
why should the second law be an exception? 



The first comment seems persuasive; however, take into account that not so academic 
models, as is the case of phantom models with constant w, do not approach equilibrium 
as they near the singularity {9]. Further, it is rather problematic to draw a dividing line 
between "academic" and "non-academic" models in cosmology. Regarding the second one, 
as far as we know, every macroscopic portion of the Universe fulfills the second law of ther- 
modynamics. Therefore, there is no compelling argument by which realistic cosmological 
models should not fulfill it as well. As for the third comment, one should not forget that 



we are dealing not with the singularities themselves but with the behavior of the models 
as they approach their respective singularities. During the approach physical laws still hold. 

From this we may conclude that the three models considered here are unphysical, at least in 
the regime nearing the singularity. The third one, the "little rip" , more particularly because 
it violates the generalized second law and it should be ruled out. 

Finally, the fact that these three models do not approach equilibrium in the last stage of 
their evolution may be seen as a feature characterizing future singularities. (This is shared 
by the phantom model mentioned above). Nevertheless, since -for the moment- the existence 
of counter examples cannot be ruled out it would be rather premature to assert the general 
validity of the said feature in FLRW cosmologies. 
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